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NLO distributions for Higgs production at the LHC 

J. Smith 3 !] 

a C.N. Yang Institute for Theoretical Physics, SUNY at Stony Brook, NY 11794-3840 USA 

We report on results for the NLO corrected differential distributions da/dpr and da/dy for the process p + p — > 
H +' X' , where pr and y are the transverse momentum and rapidity of the Higgs-boson H respectively and X 
denotes the inclusive hadronic state. All QCD partonic subprocesses have been included. The computation is 
carried out in the limit that the top-quark mass m t — * oo. Our calculations reveal that the dominant subprocess 
is given by g + g — > H +' X' but the reaction g + q(q) — > H +' X' is not negligible. Also the if -factor representing 
the ratio between the next-to-leading order and leading order differential distributions varies from 1.4 to 1.7 
depending on the kinematic region and choice of parton densities. 



1. Introduction 

The Higgs boson is the only particle in the stan- 
dard model which has not yet been discovered. 
If the Higgs mass m is between 110 GeV/c 2 and 
200 GeV/c 2 then the dominant production mech- 
anism at the LHC is g + g — > H +' X' where 
the Higgs boson couples to the gluons via a top- 
quark loop. The leading order (LO) processes 
given by g + g^ g + H, g + q(q) -> q(q) + H and 
q + q — > g + H were originally studied in [ Q , [ 
U and [ [| from which one can derive the trans- 
verse momentum (j>t) and rapidity (y) distribu- 
tions of the Higgs boson. Fortunately the calcu- 
lations simplify if one takes the large top-quark 
mass limit m t — > oo. In this case the triangle 
graphs are obtained from an effective Lagrangian 
describing the direct coupling of the Higgs boson 
to the gluons, namely 



£ eff = G<Z>(x)0(x), 
where 

0{x) = ~G^ u {x)G a ^(x). 



(1) 



(2) 



Here represents the Higgs field and G is an 
effective coupling constant which is related to the 
Fermi coupling constant Gf by 



G 2 
4^2 
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G f t 2 F 2 (t) 
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4to 2 



, (3) 



where a s {^ 2 ) is the running coupling constant 
which depends on the renormalization scale fi r . 
The function 

F(t) = 1 + (1 - t) arcsin 2 —= , (4) 
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tends to 2/(3t) in the limit of large r. Further C is 
the coefficient function which originates from the 
QCD corrections to the top-quark triangle graph 
describing the process H — > g + g in the limit 
m t — > oo. The lowest order contribution to C is 
available in [ [|] and [ || . 

A comparison was made for the differential dis- 
tributions of the LO processes in [ || where it 
was shown that the large top-quark mass ap- 
proximation is valid as long as m and pr are 
smaller than m t . The NLO matrix elements for 
g + g—tg + g + H etc., using the effective La- 
grangian were computed in [ [7J , [ || albeit in four 
dimensions. The one-loop virtual corrections to 
the LO subprocesses were presented in [g], where 
the computation of the loop integrals was per- 
formed in n-dimensions but the matrix elements 
were still presented in four dimensions. 

2. Method of calculation 

In this paper we present results from [ |l0) where 
the NLO corrections to the double differential dis- 
tributions d 2 a/dpT/dy for Higgs boson produc- 
tion in hadron-hadron collisions were calculated. 




p T (GeV) 

Figure 1. The differential cross section d a/dpT 
integrated over the whole rapidity range with 
to = 120 GeV/c 2 and /i 2 = m 2 + p\. The 
LO plots are presented for the subprocesses gg 
(long-dashed line), q(q)g (dot-dashed line) and 
100 x (qq) (dotted line) using the parton density 
set MRST98(lo05a.dat). 



Here we included all partonic subprocesses and 
used the g-g-H coupling in Eq. (3). A similar cal- 
culation has been performed in [ [II] but our ap- 
proach differs from it in various aspects. First our 
calculation is purely analytical and follows the 
calculation carried out for the Drell-Yan process 
describing vector boson production in hadron- 
hadron collisions (see [ [j^]). The approach in [ 
[j"l|| used helicity amplitudes and was based on the 
methods explained in [ [l| . Also it was mainly 
numerical and has the advantage that it can pro- 
vide exclusive distributions. In our calculation 
the matrix elements as well as the loop integrals 
and phase space integrals are computed in n di- 
mensions. Hence we could use results from [ 14 , 
[ HI, [ PI and [ O. The advantage of this ana- 



lytical approach is that one gets more insight into 
the structure of the radiative corrections. This is 
particularly important for the large corrections, 
due to soft gluon radiation and collinear fermion 
pair production, which arise near the boundary 
of phase space, where the pr of the Higgs bo- 
son gets large. Resummation of this type of cor- 
rections has been carried out for the total cross 



section in [ [lj] . Resummation of small px contri- 
butions due to the Sudakov effect has been done 
in [ [is[| . In view of the experimental problems 
to observe the Higgs boson, a recalculation of all 
the NLO corrections is necessary to be sure that 
the theoretical predictions are correct. Finally we 
mention that another paper has appeared on the 
NLO corrections to the g + g — > H + g channel, 
using the helicity framework [ ^0) . 

For our computations the number of light 
flavours is taken to be rif = 5 which holds for the 
running coupling, the partonic cross sections and 
the number of quark flavour densities. Further 
we have chosen for our plots the parton densities 
obtained from the sets MRST98 [ Jl CTEQ4 [ 
H, GRV98 [ H and MRST99 [ |J For sim- 
plicity the factorization scale fi is set equal to 
the renormalization scale /i r . For our plots we 
take p? = to 2 + p T . Here we want to empha- 
size that the magnitudes of the cross sections are 
extremely sensitive to the choice of the renormal- 
ization scale because the effective g-g-H coupling 
constant G ~ a s (fj, r ), which implies that der LO ~ 
and dcr NL ~ . However the slopes of the 
differential distributions are less sensitive to the 
scale choice if they are only plotted over a limited 
range. For the computation of the g-g-H effective 
coupling constant in Eq. (3) we chose the top- 
quark mass TO f = 173.4 GeV/c 2 and the Fermi 



constant G F = 1.16639 GeV" 



3. Results 



4541.68 pb. 



Here we will only give results for Higgs boson 
production in proton-proton collisions at the LHC 
center of mass energy y/S = 14 TeV. Since the 
hadrons H± and H2 are now identical the y differ- 
ential cross sections are symmetric (for results see 
[[hJ). In order to compare with the results in [ |ll|| 
we present LO and NLO differential cross sections 
in pt, integrated over y, for m = 120 GeV/c 2 



and pr 



p T in Figs. 1 and 2 respectively. 



The MRST98 parton densities [ |2J were used for 
these plots. We note that the NLO results from 
the qq and qq channels are negative at small pr so 
we have plotted their absolute values multiplied 
by 100. It is clear that the gg subprocess domi- 
nates but the g(g)g-subprocess is also important, 
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Figure 2. Same as Fig. 1 in NLO except for 
100*abs(gg) (dotted line) and the additional sub- 
process lOOxabs(gg) (short-dashed line) using the 
parton density set MRST98(ft08a.dat). 



yielding one-third of the total at large px- The 
LO(sum) and NLO (sum) are slightly lower when 
mt is taken very large in which case we agree with 
the results in [ |ll[ . 

There are several uncertainties which affect the 
predictive power of the theoretical cross sections. 
The first one concerns scale dependence. In the 
case of the pr-distribution one observes a small 
reduction in the scale dependence while going 
from LO to NLO. This reduction becomes more 
visible when we plot the quantity 



N 



A* 

Pt, — 
Mo 



da(p T ,n)/dp T 
da(p T ^o)/dp T 



(5) 



in the range 0.1 < /i/Mo < 10 a ^ fixed values of 
p T = 30, 70 and 100 GeV/c, see Fig.3. The 
upper set of curves at small /Vmo are for LO and 
the lower set are for NLO. Notice that the NLO 
plots at 70 and 100 are extremely close to each 
other and it is hard to distinguish between them. 
Further one sees that the slopes of the LO curves 
are larger that the slopes of the NLO curves. This 
is an indication that there is better stability in 
NLO, which was expected. However there is no 
sign of a flattening or an optimum in cither of 
these curves which implies that one will have to 
calculate the differential cross sections in NNLO 
to find a better stability under scale variations. 



Figure 3. The quantity N(pr, m/Mo) plotted in 
the range 0.1 < /J,/ Ho < 10 at fixed values of 
Pt with m — 120 GeV/c 2 and /Xq = m 2 + p T 
using the MRST98 parton density sets. The re- 
sults are shown for pt — 30 GeV/c (solid line), 
Pt = 70 GeV/c (dashed line), p T = 100 GeV/c 
(dot-dashed line). The upper three curves on 
the left hand side are the LO results whereas the 
lower three curves refer to NLO. 



The second uncertainty concerns the rate of 
convergence of the perturbation series which 
is indicated by the if-factor defined by K = 
d cr NLO /(i cr LO and finally there is the depen- 
dence of the cross section on the specific choice 
of parton densities, which can be expressed by 
the factors i? CTE Q = d a GT ^/d <r MRST and 
R GRV = d a Gm /d a MRST . The above quan- 
tities were studied in our paper, where it was 
shown that there is still a large uncertainty in 
our predictions because the if-factor varies from 
approximately 1.4 to 1.7 depending on the par- 
ton density set. The latter uncertainty is mainly 
due to the small x behaviour of the various gluon 
densities since both the partonic cross sections 
and the gluon densities increase very steeply at 
decreasing x. Finally we have defined the soft- 
plus-virtual (S+V) approximation and computed 
the ratio R s+V = da s+v / da EXACT . This approx- 
imation is quite reasonable (see Fig. 4) provided 
PT.min > 100 GeV/c in spite of the fact that x 
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Figure 4. The ratio i? s+v for the px distributions 
using the set MRST99 with y 2 = m 2 +p^ min and 
various Higgs masses given by m = 120 GeV/c 2 
(solid line), m = 160 GeV/c 2 (dashed line) and 
m = 200 GeV/c 2 (dot-dashed line). 



is still too small to belong to the large x-region. 
This means that this approximation can be used 
to resum the large corrections due to S+V glu- 
ons in order to obtain a better estimate of the all 
order corrected cross section. 
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